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Discussion of 
“MORNING-GLORY SHAFT SPILLWAYS: 
DETERMINATION OF PRESSURE-CONTROLLED PROFILES” 


by W. E. Wagner 
(Proc. Sep. 432) 


W. E. WAGNER.*—The writer is indebted to the individuals who devoted 
so much of their time and energy in preparing their constructive discussions 
of the paper. The discussions are valuable contributions to the problem of 
properly designing the overflow section of a shaft spillway. 

The graphs, prepared by Messrs. White and McPherson and supported by 
experimental data, clearly define the region of low flow in which “scale 
effects” must be considered. It is interesting that weirs having heads as low 
as 0.10 foot (figures 17 and 19) may be used without surface tension and/or 
viscosity affecting the discharge. It should be noted that these results were 
obtained in the laboratory under ideal conditions where close control of the 
approach conditions, sharpness and cleanliness of the weir, and steadiness 
of flow were maintained. However, in general practice, where ideal condi- 
tions are difficult to attain, the minimum head for flow over a weir should be 
somewhat greater, possibly as high as 0.20 foot, to assure flow conditions in 
which “scale effects” can be safely ignored, 


The discharge formula, Q = 1.66 Ve ph, developed by Messrs. 


White and McPherson for heads greater than 0.10 foot and—> ratios less than 


0.60, is a convenient method of expressing the discharge in the range of free 
flow. Discharge coefficients computed using the above formula check within 


H 
1-1/2 per cent of the coefficients shown on Figure 9 for e ratios between 


H 
0.20 and 0.50. For as ratios greater than 0.50, the deviation increases, 


H 
reaching over 4 per cent for = = 0.60, Figure 25. This deviation leads to 


the question of the upper limit at which weir flow changes in character to flow 
through a re-entrant tube. Messrs. White and McPherson showed graphically 


H 
in Figure 18 that this change occurs at 7 = 0.60. It is difficult for the writer 


to visualize an abrupt change from weir to tube flow as indicated in Figure 18 
for the following reasons: 

1. As pointed out by Mr. Blaisdell in his discussion, the head on a weir is 
usually considered to be affected when the “tail water reaches a level above 


the weir crest exceeding the critical depth of flow d..” or when i = 0.667. 


Table 4 shows that this tail water (or boil height) is reached when ts is be- 


H. 
tween 0.45 and 0.50. When = = 0.60, i = 0.94 or the high point of the boil 


is at a level above the crest equal to 0.94 of the head. A backwater of this 


* Engr., Bureau of Reclamation, Denver, Colo. 
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height surely affects the head upstream from the weir and causes a departure 
from true weir flow at some rs ratio less than 0.60. (See Figure 26 which 


H 
shows the operation of the circular weir when 7 = 0.50). 


2. In Figure 27, some of the writer’s original data are plotted on a graph 
similar to Figure 18. This plot clearly indicates that a transition zone exists 
between tube and weir flow and supports the original conclusion that the upper 

H 
limit of weir flow is 7 = 0.45 and the lower limit of tube (or “orifice”) flow 


H 
is . = 1,00. Unfortunately, actual calibration points by other experimenters 


are not available. However, the published reports by both Camp(2) and 
Gourley (3) indicate a similar transition zone since they recommend an upper 


H 
limit of 7 = 0.40 when using their formulas. 


The writer questions whether a vortex can form in the range of flows 
where a is less than 0.60 as proposed by Messrs. White and McPherson. 


At these flows, the velocity of the upper nappe surface is comparatively high 
and well directed with a clear line of demarcation between the boil and the 
nappe, as shown in Figure 26, where se = 0.50. Any vortex having a tendency 
to form in the nappe is carried immediately through the weir. Due to the 
turbulence within the boil, it is unlikely a vortex can form and maintain itself 
in this region. However, air may become entrained where the nappe joins the 
boil and perhaps affect the weir discharge, but it is questionable whether ra- 
dial piers or baffles will prevent the entrainment of air. C. J. Posey and 

H. Hsu(9) found in their vortex studies that “with purely radial inflow the 
vortex is small and transitory; its effect on discharge is negligible.” 

Mr. Blaisdell questions the writer’s statement “that most morning-glory 
spillways are designed for free flow.” Messrs. White and McPherson in their 
discussion state, “Of the 18 morning-glory spillways for which details have 
been published(3) only 4 are designed to perform submerged, 6 for ratios of 


H 
z between 0.20 and 0.40, and 8 for ratios of _§ less than 0.20.” With over 


R 
75 per cent of the reported morning-glory spillways designed to operate in 
the free-flow range, the statement can hardly be questioned. 
As suggested by Mr. Blaisdell, the discharge coefficient, Cp, in Q = Cy Ap 


H 
V2¢Hs, for zx values between 0.4 and 2.0, is given in Figure 28. For _§ 


values greater than 1.00, the coefficient is practically constant at 0.51. This 
value checks very closely the coefficient of 0.52 commonly used for flow 
through a re-entrant tube.(10) For comparison, the analytical results by 

Hs 
R 
are also shown in Figure 28. The two results compare remarkably well, 
considering the entirely different methods employed in determining the coef- 
ficients. 

The comparison by Messrs. Elassiouty, Chanda, and Johnson of the pro- 
files measured by Camp and those by the author in Figure 23 lends support 
to the reliability of the nappe shapes. It is particularly gratifying that the 
profiles check in the vicinity of the weir spring-point, since this region is 
important and extremely difficult to measure with accuracy. 


Messrs. Elassiouty, Chanda, and Johnson expressed in terms of Co and 
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feraatan by Gourley and Camp. 
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FIGURE 26- FLOW OVER CIRCULAR WEIR 
WHEN a = 0.50 
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FIGURE 27- COMPARISON OF LEHIGH AND WAGNER DATA 
IN REGION WHERE WEIR FLOW CHANGES TO TUBE FLOW 
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Discussion of 
“DIVERSION FLOW THROUGH BUFORD DAM CONDUITS” 


By Francis F. Escoffier 
(Proc. Sep. 535) 


HAROLD TULTS.*—The graphical method, developed by Mr. Francis F. 
Escoffier, is of great help in computing the water surface profiles in steady 
non-uniform flow, particularly in the open conduits of variable cross-section. 

The described application of this method to establish the reservoir levels, 
at which the flow discharge pulsates between the open-conduit-flow and the 
pressure-conduit-flow, is very instructive. However, the writer questions 
whether the predicted range of the reservoir levels, at which the pulsations 
occur, is confirmed by the actual measurements at the erected diversion 
conduits. 

There are several other factors affecting the sensitive transition of the 
flow in the diversion conduits. Besides the simplifications, made by the 
author, the following may affect the transition of the flow: (a) The inaccuracy 
of the graphical method of determining the critical depths and the minimum 
energy heads; (b) the effect of the curved flow on the pressure distribution in 
the intake transition; and (c) the variable n-value of the friction coefficient in 
the Manning formula for the part-full conduit flow. 

(a) To establish a point of tangency, presenting the critical depth for a 
certain energy level, especially at the flat F-curves, is never accurate. 
Considering the amount of the mathematical work required for the plotting of 
the F-curves of the circular and of the rectangular sections, it seems to be 
more expedient and more accurate to compute the critical depths and the 
corresponding minimum energy levels analytically, employing the general 
equation for the critical flow 


Cy coefficient of velocity (Coriolis) 
b width of water surface 
All other symbols, if not defined, have the same definition as given by the 
author. 

This Equation is adaptable for any open conduits, whatever cross- 
sectional shape it has. For a part-full circular conduit, the critical depths 
and the energy levels are easily derived by the Eq. 2 adapted for the circular 


cross-section, 


3 
Q= 


Cy Dd, -d? 


in which the area, a, of the wetted section is: 


* Asst. Hydr. Engr., Pioneer Service & Eng. Co., Chicago, Ill. 
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4 4 D 
When d-< 5 , the last two terms in Eq. 3 are negative. 


(b) The water level in an air vent, located at or near a curved boundary, 
is not determinable by the straight-line distribution of the static pressure. 
The flow of the air through the vent into the closed conduit may continue due 
to the locally reduced pressure in the curvilinear flow. Thus, a higher reser- 
voir level than that found by the author’s graphical method, would be required 
to prime the circular conduit. Of course, the pressure deviations may be es- 
tablished by drawing of the flow nets. 

(c) The application of the author’s graphical method in determining the 
water-surface profiles in the part-full closed conduits requires a modifica- 
tion of the conventional friction coefficient, which varies from the one estab- 
lished for the full conduit. 

For many years engineers(1,2,3,4) have found that using the friction coef- 
ficients, established for the pressure conduits, in the computations of the flow 
in the part-full conduits, yield incorrect results. Charles W. Sherman, (5) 
analysing the tests by O. Wilcox, found when employing the conventional hy- 
draulic radius that the friction coefficient, established for the full-pipe flow, 
should be increased up to 135 per cent for the part-full pipe. This contra- 
dicts the still wide-spread assumption that the conveyance capacity of a 
closed conduit is maximum when partly full, (for circular conduit at 
d = 0.93D). 

This may be explained by the air friction. Although not demonstrable in 
the open conduits, the retarding effect of the air friction in the part-full 
closed conduits is obviously increased by the closeness of the conduit ceiling. 

To take care of the air friction over the water surface in the part-full cir- 
cular conduit, Prof. A. Schocklitsch(6) includes the width, b, of the free water 
surface to the wetted perimeter, p, after reducing it to an equivalent wetted 
width, ab, of the same roughness as the actually wetted conduit wall. Thus, 
the total wetted perimeter of the part-full conduit is 


Popp =P EQ. 4 


Employing the same friction coefficient n (Manning) as for full conduit, the 
ratio between the mean velocities at part-full and full-pipe flow is: 


1.488 D)2/3 


2/3 
Veull 1.486 .0.5,D,2 


The quantity of { is expressed as follows (for symbols see Fig. 1): 


2/3 


In Fig. 2, the curves of the velocity ratio, Poa a 4 
Full 
coefficient,a@ , equal to 2, 1, 0.4, 0, and -0.4 are plotted. 


, for the reduction 
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The curve, a = 0, represents the variation of the velocity ratio at different 
filling depths when the same friction coefficient as for the full conduit is also 
applied for the part-full conduit, as practiced currently. 

Prof. Schocklitsch, having investigated the discharge capacities of various 
part-full circular conduits, recommends that at normal conditions when no 
forced air movements in the conduit exist, a coefficient @ equal to 0.4 should 
be used. 

When, as a result of external atmospheric conditions, an air-movement is 
induced in the conduit, the coefficient may be greater or smaller than 0.4, 
depending onthe direction of the air flow. a equal to zero is justified, when 
the air movement is in the flow direction and of the same velocity as the 
flowing water. 

As shown, the application of the author’s graphical method, excellent in 
open-conduit-flow, meets extra complications in determining the range of the 
pulsations at the transition from the open to the pressure flow and reverse, 
and also in establishing the water-surface profiles at part-full flow. 

In the solution of such delicate problems, like the one presented by the 
author, a combination of his graphical method with analytical computations 
may be helpful to increase the accuracy of the result. 
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Discussion of 
“THE LOG-PROBABILITY LAW AND ITS 
ENGINEERING APPLICATIONS” 


by Ven Te Chow 
(Proc. Sep. 536) 


SHIGEHISA IWAI.*—The paper by Prof. Chow should be highly evaluated 
in presenting quite an interesting procedure to revise Hazen’s Method with a | 
more sound theory obtained after re-examining its basic distribution which 
originally follows the log-probability law. Here, as one who is interested in 
such a research, the writer wishes to state his opinion about it as follows. 

In order to apply the log-probability law widely to the hydrologic data 
which fails in satisfying Eq. 16, we generally have to increase the number of 
parameters in addition to the two contained in the basic Equations 5 and 7. 

Prof. Chow happened to perform this treatment by adopting, not the computed 
value of Cs, but its graphically determined value. 

For the sake of such a treatment, we can, in Eq. 7, replace the reduced 
variate y = loge x by a new one represented by any of the following three 
forms: 


10g, (x + b) (i) 


y= seg, 


or as a simplified form of the above in which b = 0, 


x 
y= 106. (iii) 


where (-b) or g is a constant corresponding to the lower or the upper limit of 
the variate x, respectively, in the original frequency function (x). There- 
fore, each of Eqs. (i), (ii) and (iii) represents the distribution in which the 
value of x is restricted within the range from (-b) to (+00), or from (-b) to g, 
or from 0 to g, respectively. Also it is noted that Eqs. (i) and (iii) have one 
more additiong| parameter, while Eq. (ii) has two. 

When we take the new y by Eq. (i) and interpret Y, 7, and on the basis 
of this y, equations corresponding to Eqs. 11 and 14 are obtained: 


(iv) 
— 


* Prof. of Civ. Eng., Dept. of Civ. Eng., Kyoto Univ., Kyoto, Japan. | 
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where F 


b 
(vi) 


can be called as “a correction factor for the coefficient of variation Cy” 
similar to the form of F by Eq. 43 for C,.** 

Furthermore, it is quite interesting that Eq. 15 for C, does not change its 
form in this case and, therefore, we can obtain the following relation between 
Cg, Cy and F' instead of Eq. 16. 


Cc C 3 


Meanwhile, since the frequency factor K in Eq. 28 still maintains the same 
form as represented by Eq. 29, Prof. Chow’s principle to perform straight- 
line fitting through K is also valid for the new y. By utilizing the equation, 


6, = | log, + ( (viii) 


which is given from Eq. (v), the value of oy can be readily calculated for a 


fixed value of Cy, and for any of the assumed values of either or F', Thus, 
the obtained Gy enables us to compute K with various Ky corresponding to 
any probability P. By comparison with Professor Chow’s procedure which 
requires a rather tedious calculation of ay by Cs according to Eq. 34, the 
afore-mentioned procedure is quite simple in preparing similar tables to 
Table 2 and 4, and eventually, in constructing the chart and the probability 
paper being alike to Fig. 2 and 3. Finally, in this case, all of the values in- 


deved by Cg in these figures should be substituted by the corresponding 


——or F', while the fitting procedure remains quite the same as the author’s. 

The factor of F' can be utilized as an index to measure quantitatively and 
qualitatively the discrepancies between any estimated distribution and the 
standard one by Eqs. 5 and 7. This standard distribution should be shown as 
a straight line with an index of F' = 1 on Fig. 2. As the author stated, the 
accuracy of computed value of C, depends greatly upon the number of items 
in the series used for computing it. In this sense, the computation of Cg by 
Eq. 39, proposed again in step 1 by the author, would not be essentially need- 
ed. However, should we want to compare the computed C, with its theoreti- 
cal value, Eq. (vii) can easily give the value. 

For the distribution by Eq. (i), a solution by the method of moments was 
presented by Slade, (1’) Also, a procedure to fit a straight line, after semi- 
graphical trial caluclations, on the basis of y, was already proposed by 
Gibrat.(40) While Prof. Chow’s results in this paper would be much more 
remarkable in performing the straight line fitting on the basis of the original 


** However, for an actual adjustment enabling us to hold the relation of Eq. 


16, as shown by Eq. (vii), the reciprocal Gh) should be used to multiply 


the Cy in Eq. (v), that is computed in step 1 proposed by Hazen and Prof. 
Chow. 
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variate x, by devising a special probability paper, it should be noted that the 
above criticism by the writer is just limited on author’s process before he 
obtained the results. The laborious effort by the author in preparing the 
tables and the figures, mentioned above, must be highly appreciated. Since 
the F' and the C, in Eq. (vii) hold a one to one correspondence at a fixed Cy, 
these tables and figures will never lose their practical values. 

For the distribution by Eq. (ii), and eventually by Eq. (iii), the general 
formula of the r-th moment M, was accurately obtained by the writer, (42) 
which would be probably used for such a straight line fitting. For the distri- 
bution, especiall re Eq. (iii), the writer has developed several practical 
procedures, (2’)(3’) (4') 

By increasing the number of parameters, we can automatically improve 
the flexibility of the curve-fitting. Therefore, it is no wonder that C, of the 
extreme-value law having two parameters, is fixed at a value of 1.139 under 
varying Cy; but distributions having more than two parameters, can produce 
various combinations of Cg and Cy. However, regarding the practical flexi- 
bility of the extreme-value law, the writer has the same opinion as 
Gumbel.(56) Generally too many parameters cause a more difficult theoreti- 
cal solution and also, when we determine them by graphical inspection, too 
many combinations rather harm our quick decision on it. 

According to our experience in Japan, the extreme-value law can be ap- 
plied in a considerably wide range of Cg near Cg, = 1; and the law by Eq, (i) 
is applicable in an even wider range of C,. For the distribution extremely 
skewed, Kadoya,(5’) one of our research members, has proposed the “log- 
extreme value law” of which the basic equations are 


(ix) 


= [log ( x+b)- log, ( x,+ bd (x) 


where é is a reduced variate, Xo is the mode of x, and b, k are parameters, 
For practical usage, he has already tabulated the values of these parameters. 
It is quite interesting that this law was presented almost simultaneously with, 
but quite independently of, Gumbel’s new proposal. (6’) 

Practically, the author’s procedure in this paper would be regarded almost 
the same as the conventional methods, requiring the troublesome computation 
of moments, though it might be of great value in the United States, where 
Hazen’s Method is prevailing. We can, alternatively, take account of the 
method of maximum likelihood for the distribution by Eq. (i). In order to 
determine the three parameters contained in this distribution, a solution of 
complicated transcendental equations should be performed. However, 
Nishihara,(7’) another member of our group, has obtained very useful formu- 
las, as a result of the solution after linearizing the equations. 

The writer wishes to express his highest respect for Prof. Chow’s results 
obtained in his paper with many thanks for his kind correspondence inviting 
the writer to take part in this discussion. It is truly hoped that the further 
exchange of information in this field will be carried out between the research- 
ers in various countries. 
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R. D. GOODRICH,** M. ASCE.—For more than twenty-five years the 
writer has found that the logarithmic transformation has been “a convenient 
tool for recording, fitting, and frequency analysis of various engineering 
data.” Hence this paper has been studied with great interest. It serves to 
fill a yawning gap by presenting so ably the mathematical foundations and the 
logical developments which will make the engineering applications of the log- 
probability law more numerous and more useful. 

It is understood of course that the arithmetic mean is the most probable 
and basic value when statistical data is analyzed and assumed to be distribut- 
ed according to the normal law. When the logarithms are substituted for the 
quantities instead, their mean is the log. of the geometric mean of the quanti- 
ties making up the original data. For the benefit of the large number of engi- 
neers whose mathematical training does not extend to the heights demon- 
strated by the author, of which the writer is one, but who might wish to make 
more intelligent use of the material and methods given in this paper, it is 
suggested that the statement of the “central limit law” be amplified and more 
fully explained. Many of the references from which this information might be 
secured are not readily available except in a large university or other re- 
search library. 


H. ALDEN FOSTER,*** M. ASCE.—The writer has been interested in the 
use of log-probability plotting of engineering data for many years, going back 
to the introduction of this method for study of Flood Flows by Alien Hazen in 
1930.(39) In a discussion of Lane and Lei’s paper (59) in 1950, the writer 
pointed out the importance of having some method to determine the 


** Chief Engineer, Upper Colorado River Commission, Grand Junction, Colo. 
*** Principal Associate, Parsons, Brinckerhoff, Hall & Macdonald, New York, 
N. Y. 
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median (M) when the mean value of x (X) was known, or vice-versa. Not hav- 
ing any mathematical expression for this ratio as a function of ly, he com- 
puted the mean by arithmetic integration of the frequency function of x, -Eq. 7. 
The values of x used in this calculation were read directly from the straight 
line plotted on the log-probability paper. 

In 1952 the writer carried out additional studies of the log-probability 
function (51) in which the values of X/M, Cy and C, were determined for 
various values of Iy, again using an arithmetic-integration procedure since 
no mathematical formulas for these characteristics were available. Subse- 
quent to this study, the author’s mathematical analysis of the log-probability 
function was published.(56) When the writer compared his results for %/M 
with the values computed by the author’s formula, it appeared that there was 
a considerable discrepancy in the results. A review of the mathematical 
analysis of (56) indicated an error in the formulas derived therein, as noted 
by the author in the footnote to Eq. 10. 

The writer has reviewed the mathematical development of the present 
paper, and has verified the formulas up to Eq. 33, inclusive. As a further 
check of these formulas, he has made an arithmetic integration of the fre- 
quency curve for x, in order to obtain the characteristics X/M, C, and C,, 
using values of x taken from a table of the probability integral. On account 
of the great magnitude of the x-values at very low probabilities, points were 
selected at close spacing for probabilities less than 1.0%, the minimum spac- 
ing being 0.002%. For the rest of the curve, the points were selected at the 
mid-points of each 2% on the “%-of-time” scale. The calculations were car- 
ried out for values of ly equal to 1.00, 0.50 and 0.25. The results are sum- 
marized in Table A, in which are also shown the precise values of the same 
characteristics calculated by the author’s formulas,~—Equations 21, 14 and 
16, respectively. The latter calculations were carried out for assumed values 
of Iy, and show good agreement with the corresponding results shown on the 
author’s Fig. 1. 

Examination of Table A shows that the arithmetic integration agrees close- 
ly with the mathematical formulas for %/M and Cy except for extreme values 
of ly. The arithmetic value for Cg is also in close agreement, for ly = 0.25; 
but for larger values of Iy there are appreciable discrepancies in the Cg re- 
sults, due to the controlling effect of the extreme values of x. The writer be- 
lieves that these calculations give a satisfactory empirical verification of the 
mathematical formulas. 


Correction Factors for Computed Coefficients: 
In the writer’s paper on “Theoretical Frequency Curves” published in 
1924,* an attempt was made to establish certain correction factors to be ap- 
plied to the values of the Coefficient of Skew obtained from records of a 
limited number of items. The factors were obtained by plotting duration 
curves for certain assumed values of Cy and Cg, in accordance with the 
Pearson formula, Type II. The curves were plotted on arithmetic-probability 
paper; values of the ordinates were scaled from the diagram, for assumed 
sequences of 100, 20, 10 and 5 items; and the corresponding values of the Cy 
and Cg were computed. It was found that the Cy did not vary greatly with the 
length of the record, when the Cg did not exceed 0.60. On the other hand, 
there was an appreciable difference in the C, with records of different lengths, 
and there would be a considerable error in the C, even with a record of as 
many as 100 items. 


* Trans. ASCE, Vol. 87, 1924, Pg. 142. 
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In a discussion of the paper referred to, Allen Hazen made the suggestion 
that the value of C, computed from a record of n items could be multiplied 
by the factor, F = 1 + 8.5/n. Mr. Hazen made the same suggestion in his 
book “Flood Flows”(39) published in 1930; and this formula seems to have 
been rather extensively used since that time. The writer has felt for some 
time that Hazen’s formula was not adequate, as it was evident that the error 
in the Cg would increase for large values of Cg as well as for smaller values 
of n. The log-probability method of plotting furnished a convenient means for 
revising the correction factors. 

A series of values of the ordinates of a straight line plotted on the log- 
probability paper was obtained by use of the table of the probability function. 
The items were uniformly distributed on the %-of-time scale. 100 items were 
obtained, for 0.5%, 1.5%, 2.5%, etc; 50 items were used, located at 1%, 3%, 
5%, etc; 20 items, at 2.5%, 7.5%, etc; and 10 items at 5%, 15%, etc. Three 
different probability curves were used, with ly values of 0.50, 0.25 and 0.05 
respectively. The results were analyzed, and plotted in Fig. A. Using the 
formulas 


Adjusted Cy = Computed Cy (1 + Fy) 
Adjusted Cg = Computed Cg (1 + Fg) 


the values of Fy and F, can be taken from the diagram, corresponding to 
values of n between 10 and 100, and for various values of the Computed Cy 
and Cg obtained from the recorded items. 

It may be seen that the corrections for Cy are less than 10% if there are 
50 or more items in the record and the computed Cy is less than about 1.0; 
but the corrections for Cg become quite appreciable for all cases except 
those where the computed Cg is quite small. 

The curves of Fig. A were obtained by an interpolation process from a 
limited number of computed values; however they are believed to give an ap- 
proximate indication of the corrections which should be applied to any series 
of items which approximate a straight line on the log-probability paper. They 
should also be suitable for all practical purposes for use with other types of 
skew-probability curves. In any case, they should give more reliable results 
than the Hazen formula. 


Probability of the Mean Value: 

One characteristic of the log-probability function has not been given by the 
author. It may be convenient to determine the %-of-time at which the mean 
(%) of the x-values occurs on the probability curve. By Eq. 21, 


2 
%/M = 


Since y = logex: x = eY. Hence, & = e’m, and the median value of x, or M, = 
y where Ym = loge x. 
Hence, x/M = « el 
Or, eVm-j) = el /2 of. 
From this, (ym - y) = 1/2 (Ym - 9)/ Oy= 1/2 
Since = 2.30261, (Eq. 22a), 
 =1.1513 ly. 
The probability integral tables give values of (Ym -j)/ay vs. F(y). 
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Table A. 


Characteristics of Log-Probability Function 
By Formula Arith.-Integration 


2.30261 


1.927 1.645 


1.179 0.626 2.10 


Table B. 


Plotting Positions for Mean Value, Zs 
$-Probability 


12.4805 
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x/™ Cy Cs x/M Cy 
1.00 2.303 14.1675 13.970 2,965~ 13.620 9.87 7065 
0.75 1.726 4.328 94.0 
260 1.361 2.5969 2.397 20.94 
250 1.151 1.9401 1.663 9.59 1,0 
oS 1.035 1.7106 1.388 6.629 
4,0 0.920 1.5283 1.156 5-010 
035 0.€05 1.3837 02986 3-743 
230 0.690 1.2694 0.782 2.625 
225 0.575 1.1802 0.626 2.1.2 
220 0.460 1.1119 0-486 1.573 
015 00345 1.0615 0-356 1.113 
210 0-230 1.0269 0.233 0.712 
205 0.115 1.0066 0.116 0.350 : 
0.00 50-0 
0.10 45 4172 
0.20 40.8966 
36.4901 
040 32.2572 
0.50 
0.70 21.0148 
0.80 17.8515 
6.90 15.0062 
1.00 


for 
Coefticient of 
Voriotion 


N 
Correction Factors 


GO 70 F 


z 


COME 


WA 


fof 


Correction Factors 
for 


CoelTicien 


70 80 90 joo WN 


F/G. A 


1.0 
I 
08 
07 
04 
04 
az 
a2 
4) 
vet 
“CANIN 
JO 20 4O 
709-18 


Hence it is possible to determine the %-of-time plotting position for ym, 
which is also the plotting position for % for a given value of Iy. The computed 
values are shown in Table B. 

In this paper the author has presented an excellent mathematical analysis 
of the logarithmico-probability function, and has outlined methods for its 
practical application to enginnering problems, which should prove valuable to 
the profession. 


General 


G. N. ALEXANDER, A.M. ASCE,} and A. KAROLY.2—For some years the 
writers have made use of the log-normal distribution in hydrologic analysis. 
The comprehensive treatment by the author of this distribution and its rela- 
tion to the extreme-value distribution is particularly valuable as the statis- 
tical texts treat these distributions rather cursorily. The practical impor- 
tance of the log-normal distribution is exemplified by the author’s extensive 
bibliography. In hydrologic investigations the log-normal distribution has 
often been applied to flood flows (Author’s refs. 36, 37, 39, etc.), although 
such data is also regarded as the legitimate field for extreme-value law. 
The log-normal distribution has been less frequently applied to the daily 
stream flow data in the form of flow-duration curves, although from the 
writers’ experience on Victorian streams it is particularly applicable. It 
has therefore been used in practical applications such as estimating divert- 
ible flows for a given channel capacity, and for establishing criteria for 
stream flow accuracy. The relevant theory for these developments touches 
closely on the author’s paper and hence will be given here, but the practical 
applications will not be discussed. The writers’ discussion follows, in gener- 
al, the sequence of development used by the author. 

In the writers’ Appendix I. new symbols are suggested for the more fre- 
quent terms, with the symbols used by the author alongside. In the text of 
the discussion the author’s symbols are used, but in the writers’ Figures the 
new symbols suggested are used. 


Characteristic Relationships 


A useful addition to the relations given by the author, which the writers 
derived from Equations 5 & 6 of Kalinske’s paper (Author’s Ref. 3) is 


Cz” 


As the ratio of arithmetic and geometric means is an important parameter, 
it has been given a new symbol Cg. 

In using commercial log-probability paper (as shown in the author’s 
Fig. 2), the parameter obtained from the ratio of the ordinates at 50% and 
15.87% (or 84.13%) is og. This is therefore regarded as a more suitable 


1. Designing Engineer, Water Resources ) State Rivers and 
) Water Supply 

2. Designing Engineer ) Commission, Victoria, 
) Australia. 
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parameter than @y for graphic representation of relationships, and is there- 
fore used in the writers’ Fig. 1. For determining the other parameters a 
series of values of Cv was assumed and Cg, 2g, etc. calculated. 


Flow Duration and Volume Distribution Curves 


The flow-duration curve of the hydrologist shows the percentage of time 
(as abscissa) that the rate of flow (as ordinate) is greater or less than a giv- 
en rate. However, for many applications we are concerned with the distribu- 
tion of the volume of flow. This volume is the product of the flow rate x (in 
cusecs say) by the time during which such a rate is maintained. 

When examining this problem the writers discovered that if the time- 
distribution w(x) was a log-normal, the volume distribution was also log- 
normal, and what is more had the same standard deviation. 

The frequency function (of the flow x) from equation 73 is - 


The probability (of flow) less than the given variate is (Eq. 32) - 


k 
ae Le dky 


where Q = 1-P and is shown graphically in the lower line of the writers’ 
Fig. 2. 

The frequency function of the volume distribution is x (x) and its proba- 
bility, Q;, of volume of flow less than the given variate can be shown to be 
given by 


~ Qo 


the subscript 1 referring to the volume distribution. This equation can be 
derived from Eq. 9, by integrating between -°o and Z instead of between -co 
and +co , although it was originally proved differently. 

Thus 


and from Eq. 11 
Ky ~ Fy 


3, References to equations numbers refer to those of the author. 
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If Q; is expressed in terms of the mean x, then 


/e yr OW) d 6y) 


shown as the upper line in the writers’ Fig. 2. 
When the probabilities Q and Q2 are equal 


and hence Ky Ky 


pee 
(z,)*G 


, which is constant for any given case. 


(NOTE: Xg = Min the author’s nomenclature. In Fig. 2.C8 has the value 


4.5). Hence the two distribution lines are parallel. 
A convenient graphical method of locating the volume distribution line is 
shown in the writers’ Fig. 2, based on the relation 


=F 2c 7° OF which is derived from the previous equation. 


It was subsequently discovered that Herdan(1)4 had utilized and extended 
this principle in studying the distribution of particles; the parameter x in 
such cases measures particle size and w(x) represents the distributions. If 
x is a diametral measure, then x2 is a measure of surface area and x° of 
volume or weight. It may be shown that all the moment distributions x" y( ~) 
are log-normal with the same standard deviation of the logs. 

Thus from Eq. 10 - 


Mr ~= 


when r = 0, 1, 2, 3,etc......, we obtain 
= 1 
Taking the ratios of the moments: 


on 2 
Me, Sy 
Me 


4. Reference numbers refer to the writers’ bibliography in Appendix II. 
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“2 Cg is constant. 


To show that the respective standard deviations are constant we have 
Zz 

/ 

1)? z G 
/ 

=e 64 64-1)? X2 Cr 


and the lines representing the moment distributions are parallel. A more 
elegant proof is no doubt available. 


Relation with Extreme-Value Law 


The author gives a very useful discussion on the relation between the ex- 
treme-value law and the log-probability law, showing in his Table 3 that, 
when Cg, = 1.139 for the log-normal, the difference between the two distribu- 
tions is practically negigible. One of the writers (G.N.A.) gave a graphical 
illustration of this similarity, (2) in the case of maximum annual values of 3- 
day rainfalls in N.S.W., Australia. As the mean value of oj9 (the standard 
deviation of their logarithms) was 0.178, and hence Cg, about 1.3 it was not 
possible to discriminate on empirical grounds between the log-normal and 
extreme-value distribution. 

It is not clear what the author means by the sentence, “It means that for a 
given value of Cv = 0.364, the log-probability law furnishes a great number 
of conditions to fit curves which may be either straight or concave upward 
or concave downward as it appears on log-probability paper.” Perhaps the 
word “law” in this context is not synonymous with “distribution.” On theoret 
ical grounds, as there are two disposable parameters in each case, the log- 
normal and extreme-value distributions appear equally flexible. 

The author seeks to increase the flexibility by removing the relation be- 
tween Cy and C, in the log-normal distribution. He gives no examples and 
as yet the writers have not tried it, but for hydrologic analysis they are 
rather skeptical of the value of curve-fitting if the third moment has to be 
used. The author says that his proposed procedure has in reality a combina- 
tion of the merits of both (i.e. Hazen and Gumbel) methods. 

The writers suggest that there is a more rational way of combining the 
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2 
Hence the ratio of the means: 
x, € 
x Z2 Xn Cr 
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log-normal and extreme-value distributions, at least for hydrologic analysis. 
Thus flood flows can be regarded as the extreme-values of daily flow and as 
daily flows are often close to the log-normal distribution, flood flows in such 
cases form an extreme-value distribution derived from a log-normal parent. 
However, the extreme-value distribution used by Gumbel and upon which 
Powell’s paper is based, is the asymptotic case and makes no assumptions 
concerning the parent distribution; but as the final form is reached so slowly 
its use in practice appears unrealistic. Dr. E. J. Williams, in a private com- 
munication, shows that if the parent distribution is log-normal, the approach 
is even slower than in the case of normal distribution. For practical pur- 
poses, development of what is normally termed small sample theory is re- 
quired, but in this context the word “small” really means finite. 


Bessel’s Correction 


There are some questions concerning the application of Bessel’s correc- 
tion, i.e. the factor af . This correction has been omitted in the determina- 
tion of Cv, Cs, etc. in the section on characteristic values, and hence has not 
entered into the useful relation Cs = 3 Cv + Cv, However, in both Eqs. (38) 
and (39) the author uses the correction I . 


Using Fisher’s k statistics,(3) we have 


61 Ox? 


or since Eq. (39) can be written 


Hence Cs 


Hazen also made an adjustment to his Cs values using the factor given by 
Eq (43). 


y= 348, 
n 


It is possible that Hazen’s empirical correction may be justified in part by 
Fisher’s extension of Bessel’s correction; however as 


> 


Hazen’s correction is generally greater. 


| 

| 
| 
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CORRECTIONS 


Page 3, line 17 from top should be variance instead of standard deviation, 
and in the following line, variances instead of standard deviations. 


Page 10, the right side of equation (38) should be under the square root. 


In conclusion the writers thank the author for this paper and others he has 
written which are assisting to bridge the gap between the statistician and the 
engineer. The need for this liaison is essential in an era of specialization to 
prevent the marginal field which should be common ground becoming a no- 
man’s land. 
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Appendix I. SUGGESTED SYMBOLS 

Symbol 

DESCRIPTION 
Usec by Suggested 
Author by 
Writers 
Medien of x, or geometric meen of x M xy 
Retio of srithmetic end zeometric meens x - Cy | 

Log y Xe | 
Log j0* *10 
Mean velve of log y Xe | 
Mesn value of log 10* X10 
Stenderd devietion of x Ox o 
Standsrd devistion of loz ,x Cy Oe 
Stenderd devietion of log - S10 
Logerithm to base 10 log 10 log 
Logerithm to bese e log , ln 
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Appendix II. LIST OF REFERENCES 


1) “Small Particle Statistics,” G. Herdan, p. 118. Elsevier Publishing Co. 


2) Discussion of J. F. Mclllwraith’s paper, “Rainfall Intensity- Frequency 
Data for N.S.W. Stations,” G. N. Alexander, Journal of the Institution of 
Engineers, Australia, Vol. 26, No. 6, 1954, pp. 119-122. 


3) “Statistical Methods for Research Workers,” R. A. Fisher. Tenth Edi- 
tion, 1948. Section 14. 
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Fig.2. Derivation 
of 


Volume Distribution Curve. 
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Discussion of 
“SIMILARITY OF DISTORTED RIVER MODELS WITH MOVABLE BED” 


By H. A. Einstein and Ning Chien 
(Proc. Sep. 566) 


T. BLENCH,! M. ASCE.—The authors have discussed river model scales 
in terms of the knowledge and outlook obtained from studying bed-load trans- 
port in flumes; the writer has discussed the same subject in Separate 667 in 
terms of the knowledge and outlook obtained from studying the self-adjustment 
of large canals and rivers.(6) Superficially the results must look very differ- 
ent, because of the large number of formulas required by the one method and 
the few required by the other. Yet, if the mathematics be disregarded so that 
essential principles can be compared, the authors and the writer appear to 
have expressed identical views on all important dynamical matters, whether 
postulated or deduced. 

The concordance of essential opinions, arrived at from quite opposite view- 
points, seems to call for some explanation of the superficial differences. The 
writer prefers to develop an explanation starting from his own viewpoint of 
regime theory.(1, 2, 3) He thinks of a model-maker as being able to impose 
only three scales on a mobile-bed model; after they have been imposed all 
other scales will adjust themselves whether the model-maker wants them to 
or not. Any three independent scales are imposeable, but for the argument 
they may be taken as (i) discharge (ii) side-factor (implying tractive force 
intensity in the side zone where there is no bed-movement), and (iii) bed- 
factor (implying Froude Number in terms of depth). Obviously the discharge 
scale can be imposed directly. The side-factor scale can be imposed indirect- 
ly by calculating directly-imposeable ones from it and imposing them. The 
bed-factor scale is imposed by using a suitable bed-material, in suitable 
quantities, with a fluid of suitable nature; but, as the authors point out, and 
the writer has also stated in Separate 667, there are circumstances under 
which only unity scale is permissible. (6) To find how to produce a required 
bed-factor with a given material requires a couple of weeks of flume experi- 
ments. 

Suppose now that the model-maker had so vast a knowledge of sediment- 
transport mechanics that he knew exactly how the bed-factor depends on all 
its physical constituents. He could then dispense with the couple of weeks of 
experiment and specify, from his knowledge, that a sand of a (i) definite 
median size, with (ii) a definite dispersion of grain sizes about the median 
and (iii) a definite relative density should be used with a certain (iv) charge 
(meaning ratio of weight per second of material to weight per second of water). 
Further, the fluid should have a definite (v) viscosity and (vi) suspended load. 
To be more exact he should also specify angularity of grains, and perhaps the 
ratio of grain size to water-depth and even the ratio of channel width to depth 
in cross-over sections; but the writer believes these factors would not be 
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particularly important. Thus, for the one item bed-factor he would have at 
least six, and perhaps nine, which, added to the discharge and side-factor 
would give eleven degrees of freedom. As the writer sees the matter, the 
authors have specified most of these constituent parts of bed-factor, so have 
had to produce equations to correspond. The actual identification of the fac- 
tors from the equations is difficult without some practice, and there seems 
no point in attempting it here. 

Whether the authors’ various equations are accurate is another matter, 
and they do not appear to claim that they are. The writer feels there are two 
points concerning them that deserve comment. The first is that the equations 
rest on a prior elimination of the effect of sides by a method he thinks is du- 
bious. The second is that, according to regime theory, equation (1) contains 
D in an incorrect manner, but m ought to be 1/4. The incorrectness arises 
from the regime theory finding (from observed behavior of channels in the 
field) that the linear dimension characterising roughness should not be D, but 
a composite of the shear stress on the sides and the bed-factor; D is implicit 
in this, but would appear in the numerator instead of the denominator if made 
explicit—a paradox that used to worry river engineers till they found it was 
associated with dune movement. The writer has commented, in his discus- 
sion of Separate 611, that extrapolation of the authors’ formulas to design of 
large channels is not satisfactory; perhaps equation (1) is a principal cause, 
but recent regime analysis of flume data(4, 5) seems to show that no existing 
information possesses the accuracy, range, or completeness to justify extra- 
polation beyond a very low limit. 
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